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ABSTRACT 

We show that the orbital variation in the UV P Cygni lines of the X-ray 
binary LMC X-4 results when X-rays photoionize nearly the entire region 
outside of the X-ray shadow of the normal star. We fit models to Goddard High 
Resolution Spectrograph (GHRS) observations of N V and C IV P Cygni line 
profiles. Analytic methods assuming a spherically symmetric wind show that 
the wind velocity law is well-fit by v oc (1 — l/r) /3 , where (3 is likely 1.4 — 1.6 
and definitely < 2.5. Escape probability models can fit the observed P Cygni 
profiles, and provide measurements of the stellar wind parameters. The fits 
determine L x /M = 2.6 ± 0.1 x 10 43 erg s -1 Mq 1 yr, where L x is the X-ray 
luminosity and M is the mass-loss rate of the star. Allowing an inhomogeneous 
wind improves the fits. IUE spectra show greater P Cygni absorption during 
the second half of the orbit than during the first. We discuss possible causes of 
this effect. 
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1. Introduction 

Theorists have made considerable progress in explaining the stationary properties 
of OB star winds, and there is now a consensus that the winds are accelerated through 
resonance line scattering of the stellar continuum (Lucy & Solomon 1970; Castor, Abbott, & 
Klein 1975). The spectral lines predicted by detailed models match observations for a wide 
variety of ionization species (Pauldrach 1987), and for stars across a range of metallicity 
and evolutionary status (Kudritzki, Pauldrach, & Puis, 1987; Pauldrach et al. 1988). 

However, it has become clear that real stellar winds are time-variable and contain 
inhomogeneities which, although expected on theoretical grounds, are difficult for the 
theory to treat properly. There is independent evidence for these inhomogeneities from the 
black absorption troughs (Lucy 1982a) and variable "discrete absorption features" seen 
in wind-formed P Cygni lines (Massa et al. 1995), in the X-ray emission from OB stars 
(Harnden et al. 1979; Seward et al. 1979; Corcoran et al. 1994), and from polarimetric 
observations (Brown et al. 1995), IR observations (Abbott, Telesco, & Wolff 1984) and 
radio observations (Bieging, Abbott, & Churchwell 1989). Theoretically, the scattering 
mechanism thought to be responsible for the acceleration of the wind should be unstable, 
leading to shocks (MacGregor, Hartmann, & Raymond 1979; Abbott 1980; Lucy & White 
1980; Owocki & Rybicki 1984). Although these instabilities have been modelled numerically, 
the models encounter nonlinearities (Owocki, Castor, & Rybicki 1988). 

The unknown nature of the wind inhomogeneities has hindered attempts to infer from 
observations such properties of stellar winds as the mass-loss rates, which are critical for 
understanding massive star evolution. The observations are, in essence, averages over the 
entire nonuniform and time-variable wind, and the unknown weighting of these averages 
prevents us from knowing the magnitude of the total mass outflow. Overcoming these 
difficulties is a major goal of current studies of hot stellar atmospheres. 
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When the OB star has an accreting compact companion, X-rays from the compact 
object can be used as a probe to derive parameters of the undisturbed wind, such as its 
mass-loss rate, expansion velocity, and radial velocity law. In the "Hatchett-McCray effect" 
(Hatchett & McCray 1977; McCray et al. 1984), the X-rays remove the ion responsible for 
a P Cygni line from a portion of the wind, resulting in orbital variations of P Cygni profiles. 
This effect has been observed in several High Mass X-ray Binaries (HMXB; Dupree et al. 
1980; Hammerschlag-Hensberge et al. 1980; van der Klis, et al. 1982; Kallman & White 
1982; Haberl, White & Kallman 1989). Analysis of the velocities over which the P Cygni 
lines vary has suggested that in some cases, the stellar wind radial velocity dependence can 
be nonmonotonic (Kaper, Hammerschlag-Hensberge, & van Loon 1993). Whether this is 
intrinsic to the OB star wind or is a result of the interaction with the compact object is not 
known. 

LMC X-4 shows a more pronounced Hatchett-McCray effect than any other X-ray 
binary. The UV P Cygni lines of N V at 1240A appear strong in low-resolution HIE spectra 
at = 0, but nearly disappear at <fi = 0.5 (van der Klis et al. 1982). The equivalent width 
of the N V absorption varies by more than 100% throughout the orbit, suggesting that 
the X-ray source may remove N V from the entire region of the wind outside of the X-ray 
shadow of the normal star. 

If the X-ray ionization of the wind is indeed this thorough, then the change in the 
P Cygni profile between two closely spaced orbital phases is largely the result of the X-ray 
shadow of the primary advancing through a thin "slice" of the stellar wind. Thus models of 
the changing P Cygni profiles should be sensitive to a region of the wind that is small in 
spatial extent. Vrtilek et al. (1997, Paper I) obtained high- resolution spectra of the N V 
lines (and the C IV AA 1548, 1550 lines, which show weaker variability) with the Goddard 
High Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope. These spectra, 
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unlike the earlier IUE spectra, provide a high enough resolution and signal to noise ratio to 
allow us to examine the detailed velocity structure of the P Cygni lines over closely spaced 
intervals of the binary orbit. 

LMC X-4 consists of a main-sequence O star (variously classified as 07-9 III-V) of 
about 15 M and a ~ 1.4 M neutron star with a spin period of 13.5 seconds (Kelley et 
al. 1983). Pulse delays from the neutron star have established the projected semimajor 
axis of the orbit precisely as 26 light seconds. The 1.4 day orbital period is also seen in 
X-ray eclipses (Li, Rappaport, & Epstein 1978) and modulations of the optical brightness 
(Heemskerk & van Paradijs 1989). The X-ray spectrum in the normal state is hard (photon 
power-index a 0.8) with little absorption of the soft X-rays (Woo et al. 1996 give 
N H = 1.14 ± 0.005 x 10 21 cur 2 ). The total quiescent X-ray luminosity is ~ 10 38 erg s -1 , 
near the Eddington limit for a 1.4 M Q neutron star. 

X-ray flares, in which the pulse-averaged X-ray luminosity can exceed 10 39 erg s -1 , 
occur with a frequency of about once a day (Kelley et al. 1983; Pietsch et al. 1985; Levine 
et al. 1991). These flares have been interpreted as dense blobs of wind material crashing 
down on the surface of the neutron star (White, Kallman, & Swank 1983; Haberl, White, & 
Kallman 1989). It is an open question whether these blobs represent the natural state of the 
stellar wind, or are produced by interactions between the X-ray source and an unperturbed 
stellar wind. An accreting compact object embedded in a stellar wind can affect the wind 
through the heating, ionization, and radiation pressure of its X-ray emission, and through 
its gravity. Numerical simulations (Blondin et al. 1990, 1991) of the modifications of the 
stellar wind by the compact object show that the wind can form accretion wakes and 
disk-like structures (even in systems which are not thought to have substantial mass-transfer 
through Roche lobe overflow). Thus an investigation of instabilities in the wind of LMC X-4 
also bears on the cause of the X-ray flares. 
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We present a series of increasingly more sophisticated models of the P Cygni lines in 
LMC X-4, as observed with the GHRS and reported in Paper I. In §^|, we compare the 
N V and C IV P Cygni lines to those of isolated OB stars and to those of other massive 
X-ray binaries. We make inferences from analytic and approximate models to the P Cygni 
line variations in §|3|. Then we attempt more detailed models in §[|, and derive best-fit 
parameters for the wind and X-ray luminosity. 

We note that this is the first detailed fit of P Cygni lines that has been attempted for 
a HMXB. The sophistication of the models that we present is justified by the high quality 
of the data on this system that should be obtained in the near future. 

2. Comparison spectra 

In this section, we compare the UV line profiles of LMC X-4 with the line profiles of 
similar stars. 

2.1. Other X-ray binaries 

The strength of the orbital variation in P Cygni lines (the Hatchett-McCray effect) 
differs among the HMXB. Figure 1 compares the N V lines at two orbital phases with the 
corresponding features in 4U1700-37 and Vela X-l (the latter two objects were observed with 
IUE at high spectral resolution in the echelle mode). The optical companion of 4U1 700-37 
is very massive (~ 50 M ) and hot (spectral type 06.5 with effective temperature 42,000K; 
Kaper et al. 1993). In this system the N V resonance line is saturated and shows little or 
no variation over orbital phase; this saturation is caused by the very strong stellar wind 
produced by the extremely massive and hot companion. The companion of Vela X-l has 
a lower mass (23 M Q ) and is considerably cooler (spectral type BO. 5 with an effective 
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temperature of 25,000K; Kaper et al. 1993). The wind in LMC X-4 is expected to be 
weaker than in Vela X-l and 4U1 700-37, as the optical companion is not a supergiant. As 
the abundances of metals are lower in the LMC than in the Galaxy, lower mass-loss rates 
for a given spectral type are expected. The low abundances will also make the N V lines 
less saturated and the effects of X-ray ionization more visible. The X-ray luminosity of 
LMC X-4 is the highest of the three systems. As a result of these considerations, it is not 
surprising that the effects of X-ray ionization on the N V lines are strongest in LMC X-4. 
For all three systems the C IV resonance lines show little orbital variation during the phases 
observed. 



2.2. Photospheric lines at ~ 0.5 

The effects of the X-ray photoionization on the blueshifted P Cygni absorption are 
greatest for LMC X-4 near <p = 0.5. The observed N V line profiles still show residual 
absorption at = 0.42. Paper I interpreted this absorption to be the photospheric N V 
doublet, which is masked by the wind absorption during X-ray eclipse. If this interpretation 
is correct, the width of the absorption is determined by the rotation of the normal star, 
not the wind expansion velocity. Alternately, at the base of the wind where the expansion 
velocity is low, the wind may be dense enough that N V can still exist even in the presence 
of X-ray illumination. 

Nearly all early-type stars whose N V photospheric absorption lines could be usefully 
compared with those of LMC X-4 also have strong wind absorption. The 09 IV star //Col 
(HD 38666) seems to be a promising candidate for comparison, as its N V line has negligible 
red-shifted emission, and the blue-shifted absorption is confined to features near the rest 
wavelength. 
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Although this is the best candidate for the intrinsic N V photospheric absorption 
spectrum that we were able to obtain, Wollaert, Lamers, & de Jager (1988) suggest that 
the N V and C IV lines in this star are wind-formed because they are asymmetric. We have 
compared IUE and HST spectra of /iCol with our spectra of LMC X-4, and although the 
detailed fit is poor (Figure 2a), the strength of the absorption in /iCol is comparable to the 
absorption in LMC X-4 at <fi = 0.4. This provides qualified support that the photospheric 
N V lines of subgiant O stars can be comparable to the lines we see in LMC X-4. 

We also compare the LMC X-4 spectra with the optical absorption lines which 
Hutchings et al. (1978) concluded were photospheric. Hutchings et al. found that the 
optical lines were consistent with an orbital velocity of 50 — 60 km s~ x and a systemic 
velocity of 280 km s" 1 . They measured the width of the optical lines to be ~ 170 km s _1 , 
which would be expected if the primary rotates with a period twice as long as the binary 
orbit. We find a good fit to the N V absorption profile using Gaussians with width 
(approximately equal to vsini width) 170 km s" 1 and optical depths Tb = 0.6 and r r = 0.3 
in the blue and red components (Figure 2b). 

While the N V absorption velocity agrees with both the optical absorption line velocity 
found by Hutchings et al. and with the systemic velocity of the LMC (Bomans et al. 
1996), the cross-over between absorption and emission occurs at relative red-shifts of 
100 - 200 km s -1 . This may suggest that LMC X-4 has a red-shift relative to the LMC and 
that the "photospheric" absorption is actually blue-shifted P Cygni absorption. 

In the absence of a clear means of separating the wind and photospheric absorption, 
we will allow, in the more detailed analysis that follows, the strength of a photospheric 
absorption line to be an adjustable free parameter. We note that Groenewegen & Lamers 
(1989) fit P Cygni lines observed with IUE, and that for the stars HD 24912, 36861, 47839, 
and 101413 (with spectral types 07.5111,08111, 07 V, and 08 V respectively), the best-fit 
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values of r r for N V are 0.0, 0.4, 0.3, and 0.3 respectively. 



3. Simple Models 

In the discussion that follows, we assume that the orbital variation of the N V 
P Cygni lines is caused principally by the Hatchett-McCray effect, that is, through the 
photoionization of the stellar wind by the X-ray source. There do not appear to be any 
variable Raman-scattered emission lines which Kaper et al. (1993) suggested cause an 
orbital variability in the X-ray binary 4U 1700-37. 



3.1. Equivalent width variations (IUE and HST observations) 

Following van der Klis et al. (1982), we examine the orbital variation of the equivalent 
width of N V in LMC X-4 using a simple numeric model. Hatchett & McCray (1977) 
showed that the X-ray ionized zones could be approximated as surfaces of constant q, where 
each point in the wind has a value of q given by 

= (1) 

Here, L x is the X-ray luminosity, n x is the electron density at the radius of the X-ray 
source (D = 1.77 stellar radii in the case of LMC X-4), and £ is the ionization parameter 
£ = L x /n c r^, where n e is the electron density at the given point, and r x is the distance from 
that point to the X-ray source. Given a constant mass-loss rate M for the wind, n e can be 
found from mass conservation: 

M -s 

ne = L2 4^; m « ■ (2) 

We illustrate these parameters in Figure 3. 

For large values of q, the surfaces of constant ionization are spheres surrounding the 
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X-ray source, while for q < 1, the surfaces are open. Where the surfaces extend into the 
primary's X-ray shadow, a conical shape must be used instead of the constant q surface. 
We model the absorption line variation using the velocity law v oc (1 — 1/x), where x is 
the wind radial distance in units of the stellar radius. According to van der Klis et al. the 
results are not sensitive to the velocity law assumed. We account for the orbital inclination 
i (i = 66.2li;|°, Pietsch et al. 1985) by using a corrected orbital phase 0' such that 
cos27T0' = cos27T0sku. To model the orbital variation of equivalent width, we used the 
best-fit escape probability model we describe in §4.1|, except that the only X-ray ionization 



occurs within the zone bounded by q, which is completely ionized. 

We measured the equivalent width of N V, C IV, and Si IV for all the IUE observations 
of LMC X-4 (using the Final Archive version of the data, which has been processed by 
NEWSIPS, the NEW Spectral Image Processing System, de la Pena et al. 1994). We 
used the ephemeris of Woo et al. (1996) to find the corresponding orbital phases for each 
spectrum. In Figure 4, we compare theoretical variations in equivalent width with those 
observed with IUE and the GHRS. That is, we plot 

A = (W\ — W x )/W x (3) 

where W\ is the equivalent width, and W\ is the average value of W\. Here, W\ measures 
only the equivalent width of the absorption, and not the P Cygni emission. As the 
GHRS spectra were obtained during less than a complete binary orbit, we use the average 
equivalent width from the IUE observations to find the per cent variation in equivalent 
width for the GHRS spectra as well. 

The results shown in Figure 4 are consistent with q ^ 1, in which case the X-ray 
source photoionizes at least the entire hemisphere opposite the primary, and possibly as 
much as the entire wind outside of the primary's shadow. The fit to the IUE data is better 
for q = 1 (x 2 = 5.5) than for full X-ray photoionization (x 2 = 11), although the GHRS 
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equivalent widths agree with the prediction of full ionization. We note that interstellar 
and photospheric absorption lines, if present, provide additional absorption components 
that are not expected to show strong orbital variation. If these components are present, 
this would naturally explain why the fractional variation in equivalent width that we see 
is less than that predicted by q <C 1. For example, for N V, W\ = 2.4A, while W\ due to 
photospheric absorption was found with the GHRS to be ~ 0.8A. Without consideration 
of the photospheric absorption, the observed variation of 1.6A corresponds to A = 67%. 
However, this corresponds to a 100% variation in the equivalent width of absorption due to 
the wind alone, excluding the photospheric absorption. 

It is likely that Figure 4 underestimates the per cent variation of the P Cygni absorption 
for the C IV and Si IV lines, as X-ray should be more effective in removing these ions from 
the wind than removing N V from the wind. A systematic uncertainty in choosing the 
continuum level will have a greater effect on the fractional variation in the equivalent width 
of the weaker lines, such as Si IV. 

In the next section, we will make the simplifying assumption that the entire region 
outside the conical shadow of the normal star is ionized by the X-ray source, and from this 
derive the radial velocity law. 

3.2. The radial velocity law from the maximum velocity of absorption 

The spectra obtained with the GHRS were read out at 0.5 second intervals, using 
its RAPID mode (Paper I). The P Cygni absorption profile was variabile on time scales 
shorter than the HST orbit. We interpret this as the result of the X-ray shadow of the 
primary crossing the cylinder subtending the face of the primary, where the blue-shifted 
absorption is produced. We attempt to infer the radial velocity law in the wind from this 
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rapid P Cygni variation and a few simple assumptions. 

We assume here that X-ray photoionization removes N V and C IV from the entire 
wind outside of the X-ray shadow of the primary. This assumption is easily shown to be 
reasonable, based not only on the IUE and HST equivalent width variations (§ |3.1| ) but also 
on the expected value of q. Given a mass-loss rate of 10 -6 M yr _1 , a velocity at the X-ray 
source of 600 km s -1 , and an ionization parameter log£ = 2.5 at which there should be 
little N V present (Kallman & McCray 1982), we find q = 0.16. 

As Buff & McCray (1974) pointed out, as an X-ray source orbits a star with a strong 
stellar wind, the absorbing column to the X-ray source should vary in a predictable way that 
depends on the wind structure and the orbital inclination. If the wind is not photoionized, 
then the absorbing column increases near eclipse. In a companion paper (Boroson et al. 
1999) we will examine the ionization in the wind from the variation in the absorbing column 
as seen with ASCA. Preliminary analysis supports q <C 1. 

We now use the geometry shown in Figure 5. The primary star is assumed to be 
spherical with a radius R. The neutron star orbits at a radius D. The angle between the 
ray connecting the neutron star to the center of the primary and a tangent to the primary's 
surface is given by cos 8 = R/D. Then we define three relevant distances: 

COS(7T — 6 — ZTTCp) 

D 3 = D 2 + R ° 2 +— g— , (5) 
D sin 2n(j) tan 2n(f) 

D 4 = (D 2 3 + R 2 ) 1/2 . (6) 

We wish to compute, given 0, the radius r(0) in the wind responsible for the maximum 
observed velocity of absorption. Given (r ((/)), v (([))) derived from closely spaced when 
the shadow of the primary is moving across the line of sight, we can compare with the 
theoretical expectation of the wind acceleration velocity law. (Note that this assumes a 



-13- 



spherically symmetric wind.) 

Given a monotonic velocity law, in most situations the maximum absorption velocity 
will be observed to be V(4>) corresponding to a wind radial velocity V milx (Figure 5), at a 
distance of D 4 from the center of the primary. This is the furthest radial distance at which 
the wind can absorb stellar continuum and at which the X-rays from the neutron star are 
shadowed by the limb of the primary. Note that the wind at this point is moving obliquely 
so that V(4>) = V max D 3 /D 4 . As increases from 0.1 to 0.5, the wind at V mSLX will move 
with increasing obliqueness to the line of sight. As a result, even though this point has 
the greatest radial distance of any point in the wind that can absorb N V, it still may not 
cause the absorption with the highest possible line of sight velocity, even in a monotonic 
wind. As a first check to see if the obliquity of the velocity vector at V max is important, we 
can compare V(<f>) with Vjine— of —sight 5 which can be calculated from D 2 , assuming a radial 
velocity law. 

We assume that D/R = 1.77, which can be calculated from the mass ratio and the 
assumption that the primary fills its Roche lobe. (Note that we use the separation between 
the centers of mass of the two stars, and not the distance of the NS from the CM of the 
system, which is about 10% smaller.) 

We determine the maximum velocity of absorption by fitting the high-velocity edge of 
the absorption profile with a broken line (Figure 6). The line is horizontal for the continuum 
blue-ward of the maximum velocity of absorption and has a nonzero slope, allowed to vary 
as a free parameter, at velocities less than the maximum. When the doublets do not overlap 
(<f) > 0.19), we are able to obtain two maximum velocities at each phase, one for each 
doublet component. For the red doublet component, we do not fit a broken line in order to 
determine the maximum absorption velocity, as the profile is never horizontal. Instead, we 
find the maximum absorption velocity in the red doublet component by determining where 
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the flux rises to the continuum level (which is determined by the broken-line fit to the blue 
doublet component.) The results from the red and blue doublet components do not show a 
systematic difference. Thus averaging maximum absorption velocities determined in each 
line allows us to reduce the errors in our velocity measurements. 

We can apply a further refinement to the method. The observed P Cygni trough 
is actually a sum of the true absorption profile and blue-shifted P Cygni emission. The 
observed maximum absorption velocity may be shifted by the addition of the emission. 
Although we can not observe the blue-shifted emission profile directly (it is masked by the 
absorption whose profile we are trying to uncover) we can approximate it using the reflected 
red-shifted profile at a complementary orbital phase, assuming spherical symmetry and the 
Hatchett-McCray effect. The emission flux F at velocity v and at phase is given by 

F(v,(j>) taF(-v, 0.5-0). (7) 

For example, at = 0.1, the red-shifted emission at high velocities will be diminished as a 
result of photoionization; the blue-shifted emission at = 0.4 will be similarly diminished. 
(This symmetry is broken because the primary shadows a portion of the red-shifted emission 
but none of the blue-shifted emission. However, we can safely neglect this effect because 
the blue-shifted emission in front of the stellar disk can not have a higher velocity than the 
blue-shifted absorption.) 

The zero-velocity for the reflection of the red-shifted emission is not completely 
constrained. We allow two choices: 1) that the systemic velocity of LMC X-4 is identical to 
that of the LMC, and 2) that the systemic velocity of LMC X-4 is 150 km s" 1 so that the 
turnover between red-shifted emission and blue-shifted absorption occurs near a velocity 
of 0. 

We give the measured values of (r(0),t>(0)) in Table 1. We consider that a reasonable 
la - error for the velocity measurements is 50 km s _1 . 
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As an independent check on our measurements of wind velocity, we examine the 
maximum velocities at which the absorption profiles vary between two orbital phases. This 
is the only method we can apply to the C IVA1548 lines, in which the P Cygni absorption 
can not be separated from the strong interstellar and photospheric lines by examination of 
the spectrum alone. The best-fit maximum velocity of the C IV line is at 910 ± 50 km s _1 . 
We note that we did not observe any red-shifted P Cygni emission in C IV, so were unable 
to compensate for the variability in the blue-shifted emission. 

The maximum velocity of the change in N V absorption between (ft = 0.111 and 
(ft = 0.209 is 1130 ± 50 km s _1 , consistent with the measurements in Table 1. Between 
(ft = 0.209 and (ft = 0.312, the peak velocity of absorption variations is 560 ± 50 km s _1 , 
consistent with the previous measurements (depending on the model of the blue-shifted 
emission that is subtracted). However, between eft = 0.312 and (ft = 0.410, the maximum 
velocity of the absorption variations is 160 km s _1 . This is lower than the velocity inferred 
from the line profiles alone (~ 320 km s _1 at (ft Pd 0.312). 

As discussed in § |2~^ , the blueshifted absorption may have a component which is 
produced near the stellar photosphere, rather than in the wind. In this case, the velocity 
width of the line may be produced by the rotation of the star, rather than the wind 
expansion. We suggest that the absorption line for (ft > 0.25 is likely to contain a large 
photospheric component. Thus we do not fit to the maximum absorption velocity at these 
phases. 

We fit our measurements of (r((ft),v((ft)) to the analytic expression 

v = v oo (l-l/x) p + v (8) 

where Vq is the systemic velocity of LMC X-4, and x = r/R. The best-fit values and formal 
errors are presented in Table 2, for 3 possibilities of organizing the data, depending on 
whether we subtract our model for the blue-shifted emission from the spectra and whether 
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the velocity of the reflection point is or 150 km s -1 . We find the best fit using a downhill 
simplex method (Nelder & Mead, 1965) and estimate errors using a Monte-Carlo bootstrap 
method. 

We plot (r ((f)), v ((f))) (measured from both doublet components) in Figure 7, along 
with the analytic fits. There is no evidence that the wind is nonmonotonic. All the results 
indicate that the wind expands more slowly than expected from the radiatively accelerated 
solution (which can be well-approximated with (5 = 0.8; Pauldrach, Puis, & Kudritzki 
1986). While there is a large range of (3 that can fit the 3 arrangements of observing 
(r (0), v ((f))) (1.4 < f3 < 2.4), the choice of (5 — 1.39 ± 0.14 provides the lowest \ 2 and 
the most reasonable systemic velocity. (However, note that turbulence in the wind could 
add a constant velocity to the absorption edge and could cause the data to be well-fit by 
a nonphysical value of the systemic velocity.) The models are also consistent with the 
measurements of the wind veocity from differences of successive spectra of C IV at = 0.15 
and N V at = 0.3. 

While a coincident stellar absorption feature at A m 1249A decreases the red-shifted 
emission for velocities v ^ 1000 km s" 1 , compensating for this feature in our subtraction of 
the blue-shifted emission would only increase slightly the value of (5 that we measure. 

Our fits for the stellar wind acceleration parameter (5 are not sensitive to our choice of 
D = 1.77; using D = 1.70 and D — 1.84 we find f3 within la of values in the first two rows 
of Table 2, but f3 = 2.00 instead of 2.37 ± 0.11 for the third row. 

The effects of microturbulence in the wind on the measured maximum wind velocity 
should not alter the values we have found for (3. For example, if the microturbulence 
velocity is constant, then our fits would give a different value of vq but the same (3. If 
the microturbulence velocity is a constant fraction of the wind velocity, then as a rough 
approximation, all measured values of v((p) will be multiplied by a constant, leading us to 
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fit (r ((f)), v ((f))) with a different value of but the same value of f3. 

In this section we have demonstrated a new method of probing the velocity field 
of a stellar wind. We have assumed only standard parameters of the stellar radius and 
separation, that the wind is spherically symmetric and that the X-ray source completely 
ionizes the wind. One drawback of the method is that only a very small portion of the 
available data is used, that is, only the velocity of the edge of the absorption. In the next 
section, we adapt the method to fit the entire absorption profile. 

3.3. The Radial Velocity Law from the Absorption Line Profile 

We now attempt to fit the entire absorption profile to determine the wind radial 
velocity law. The geometry of this method is shown in Figure 8. The X-ray shadow of the 
primary is conical. At distances r ^> R from the center of the primary, the wind within the 
cylindar subtending the stellar surface has only a small oblique velocity (assuming a radial 
velocity field). Thus the surfaces of constant projected velocity, far from the primary, are 
planes perpendicular to the line of sight. The intersection between the cone of shadow and 
the planar isovelocity surface at a given observed velocity will most often be a parabola. 
This parabola bounds those portions of the stellar disk that lie behind unilluminated gas 
at the observed velocity, and thus are absorbed by the N V in the wind, and those portions 
that do not. 

If we could obtain an image of the stellar disk of the primary in LMC X-4 with a 
narrow-band filter centered at a wavelength corresponding to a Doppler shift at the given 
velocity, we would see a crescent-shape bounded by the star's disk and the parabola. As 
the wavelength of the filter approached the terminal velocity of the wind, the intersection 
of the plane and cone would bound a smaller and smaller region of the stellar surface, and 
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would thus determine the shape of the edge of the P Cygni profile. 

In general, the isovelocity surfaces will not be planes, and we will need to solve for 
them numerically. We can no longer simply read off the velocity corresponding to each 



radius from the absorption profile, as we did in § |3.2| . Instead, we need to assume a velocity 
law, find the isovelocity surfaces, find their intersection with the conical shadow, and then 
compute the absorption profile and compare with the observed profile. 

To use the observed absorption profile, we must correct for the blue-shifted emission by 
subtracting the red-shifted emission at a complementary phase. Now that we consider the 
line profile and not merely the edge of the absorption, we must recognize the asymmetry 
between the red and blue-shifted emission profiles at <ft an d 0.5 — <ft which is introduced by 
the primary's shadow of the red-shifted emission. We correct for this by assuming that the 
shadowed red-shifted emission at the complementary phase is given by 

F shadow (v, <j)) = ~ [\ - VX2 ~ 1 ) F,U-v, 0.5 - 0) (9) 

where x is the radial distance in the wind, normalized to the stellar radius, -F shadow is the line 
profile of the shadowed part of the red-shifted emission, and is the intrinsic absorption 
line profile. For this approximation, we have noted that the shadowed red-shifted emission 
at 4> corresponds to the blue-shifted absorption at 0.5 — <p, with the dilution factor needed 
to compare the absorption with the emission. 

We also need to assume an optical depth in the wind. As a first attempt, we assume 
that the wind is entirely black to the continuum in the region of N V. We neglect limb 
and gravity darkening of the primary, and its nonspherical shape. We correct for the 
photospheric absorption by subtracting the profile at <ft = 0.4, when we assume the 
absorption is entirely photospheric. 



A fit to individual short subexposures, as in §3.2 would not be time-efficient. Instead 
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we fit to the blue doublet component absorption profiles at = 0.1, 0.2, 0.3. The fit is poor, 
with xl — 11 with 253 degrees of freedom, = 1780 ± 10 km s -1 , v = 110 ± 10 km s _1 , 
and f3 = 4.25 ± 0.08, which we do not consider a reasonable value. Allowing a uniform 
optical depth r in the wind gives a better fit, with x 2 — 2.95, Voq = 1330 ± 30 km s _1 , 
v = 72 ± 24 km s" 1 , (3 = 1.57 ± 0.04, and r = 0.97 ± 0.06. We show the results of this fit 
in Figure 9. This high value of (3 is similar to those found in §|3T2". 



In the following section, we will investigate detailed numeric models of the wind, 
relaxing our assumption that the X-rays ionize the entire wind outside the X-ray shadow of 
the normal star. We do not remove the restriction that the wind velocities are spherically 
symmetric, as this would allow more than one point along the line of sight to be at the same 
projected velocity, which would require us to use more sophisticated but time-consuming 
techniques, such as Monte-Carlo simulation of P Cygni profiles. 



4. Escape Probability Models 

4.1. Methods and basic models 

To model the P Cygni line variations, we apply and extend the method of McCray et 
al. (1984). We compute the line profiles using the Sobolev approximation and the escape 
probability method (Castor 1970; Castor & Lamers 1979), using a wind velocity law given 
by Equation || with [3 = 1. The optical depth in the undisturbed wind obeys 

r = T(l + 7 )r" 7 (10) 

where T and 7 are empirical parameters (Castor & Lamers 1979). 

The optical depth of gas illuminated by the X-ray source is given by the Sobolev 
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approximation, 

t = (ne 2 1 mc) f \ riij(dv / dr)^ 1 , (11) 

where e is the electron charge, m the electron mass, Ao the rest wavelength of the transition, 
/ the oscillator strength, and n^j the particle density of element i in ionization state j. 
From mass conservation, n it j is given by 

M 

niJ = 4^ mR a ' l%h (12) 

where gij is the fraction of element % in ionization state j, and is the fractional abundance 
of element i. Given a value of £ and the shape of the X-ray spectrum, the ion fraction 
g iy j (and thus ny) can be computed. For this we used the XSTAR code (Kallman & 
Krolik, based on Kallman & McCray 1982), which requires that the X-ray luminosity L x 
be integrated over the energy range 1-1000 Ry (13.6 eV-13.6 keV). The effectiveness of 
X-ray photoionization depends on the X-ray spectrum; we have used a broken power-law 
spectrum, 

F(E) = KE - 37 for E < 25 keV (13) 
= KE~ 2 for E > 25 keV, (14) 

where F(E) is the energy spectrum. The power-law index for E < 25 keV was determined 
from simultaneous ASCA observations of the 0.4-10 keV spectrum (Paper I), while the 
high-energy cutoff has been observed at other times using GRANAT (Sunyaev et al. 1991). 
Adding a soft blackbody that contributes ~ 10% to the total flux (as seen in the X-ray 
spectrum with ASCA) has little effect on the ion fractions. We can treat the combined 
effect of the background ionization (due to the photospheric EUV flux or X-ray emission 
from shocks in the wind for example) and the X-ray ionization by 

SUtot = (l/sux + VsUback - I)" 1 , (15) 
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in the case in which the wind in the absence of X-ray illumination is at least as highly 
ionized as the ion which causes the P Cygni line. If this is not the case, then we have 

l/fl't.i.tot - 1 = (l/#i,j~l,X + V&J-l.back - 2) _1 (16) 

In Equations [15] and [16|, g^x and flyback are ; respectively, the ion fractions with X-ray 



illumination (computed using XSTAR) and without (computed from equations [10|, [TT], and 



We assume that all N in the unilluminated wind is ionized to at least N V and apply 



Equation yj. The models of rSco (T ^30,000) by MacFarlane, Cohen, & Wang (1994) 
predict that the dominant ionization stage of N in the wind is N VI for wind velocities 
v > 0. 21*00. In these models, the wind is ionized by X-rays, presumably caused by shocks in 
the wind, that have been observed by ROSAT. However, MacFarlane et al. found that their 
model of rSco did not match UV observations, and that X-rays were less effective in ionizing 
the denser winds of hotter stars (consistent with the results of Pauldrach et al. 1994). We 
also tested models in which N IV was the dominant ionization stage of the unilluminated 
wind, and found that the best-fit parameter values were within the errors of the model we 



present here. The choice of Equation [15] or Equation [16] has a small effect on our results, 
as it affects only that portion of the wind that is illuminated by direct X-rays from the 
neutron star, yet is ionized so slightly that the background ionization rate is important. 

We treat the overlap of doublet components using the approximation of Castor & 
Lamers (1979). The optical depths in the photospheric absorption lines are assumed to 
have Gaussian profiles. The red and blue doublet components of the photospheric lines have 
optical depths at line center of Tb and r r = 0.5Tb ( m the lines we consider, the oscillator 
strengths are in a 2:1 ratio). 

We do not account for limb and gravity darkening of the companion or its nonspherical 
shape. However, we normalize the absorbed continuum (which arises in the visible face of 
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the star) to the continuum level observed by the GHRS at each particular orbital phase, 
while we normalize the scattered emission (which can originate anywhere on the companion) 
to the average continuum level observed with the GHRS. 

The free parameters of the fit are T, 7, v^, v (the systemic Doppler shift of LMC X-4), 
L 38 /M_ 6 (the ratio of X-ray luminosity in 10 38 erg s _1 to mass-loss rate in 10~ 6 M yr -1 ), 
M_ 6 , and r b . We employ a downhill simplex method (Nelder & Mead, 1965) to find the 
minimum x 2 and then refine the solution and find errors in the parameters using the 
Levenberg-Marquart method (Marquardt, 1963). For all velocities, we set the error to at 
least 35 km s" 1 , which is the resolution limit of the GHRS G160M grating. We exclude 
from the fit the central lA region around the rest wavelength of each transition, which is 
poorly modeled by the escape probability method. 

In Table 3 we show the best-fit parameters and we display the fits to the line profiles 
in Figure 10. 

The fits determine L 38 /M_ 6 = 0.26 ± 0.01 because the ion fractions in the wind 
(and thus the optical depths) depend sensitively on this ratio. The fits also determine 
M—qCLn /on,lmc — 3.2 ± 2.3, where on.lmc — 1-8 x 10~ 5 is the Nitrogen abundance in the 
LMC, assumed to be 20% of the cosmic value. The error determined for this parameter 
shows that it is not effectively constrained by the fit; this is because it only affects linearly 
the optical depth in the portion of the wind outside of the X-ray shadow. 

In general, the model successfully matchess the change in the P Cygni profiles between 
successive orbital phases (Figure 10, right-hand panel). However, the predicted change in 
the N V profiles between = 0.31 and <fi = 0.41 is blue-shifted from the observed change. 
This is further evidence that the wind expands more slowly at small radii than expected 
from a (5 — 1 velocity law. 
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4.2. Wind Inhomogeneities 

We now allow the wind density to be inhomogeneous by adding two new free 
parameters, another mass-loss rate ZM, where Z > 1, and a covering fraction /. In the 
absence of X-ray photoionization, the optical depth in the two wind components are t% — r 
(where r is given by Equation [10]) and T2 = Zt. The ionization fractions are computed 
for a given point in the wind that has a density determined from M and ZM. For the 
blue-shifted absorption, a fraction / of sightlines to the stellar surface encounters material 
with the higher density. This material has a much higher optical depth, as not only is there 
more of this material, but it also contains a higher fraction of N V. 

Applying a downhill simplex x 2 minimization algorithm gives an improved xl = 4.01. 
The fit is improved for two reasons. First, the fit of the homogeneous wind model 
over-predicted the effects of X-ray ionization on the red-shifted emission (Figure 10). The 
inhomogenous wind model allows more N V to survive in the X-ray illuminated wind, 
so that more of the red-shifted emission persists. Second, the homogeneous wind model 
predicts too much absorption in the 1238A line and not enough absorption in the 1242A line 
at <p = 0.12. When an inhomogeneous wind is allowed, the strength of the absorption in the 
two components can be more nearly equal, as absorption from the dense wind is saturated. 
We compare the inhomogeneous wind models with the GHRS observations in Figure 11. 

We find a "covering fraction" / = 0.33 ± 0.01; this is the probability that a ray from 
the primary encounters denser material with the mass-loss rate ZM. Our best-fit values 
are Z = 700 ± 200 and Z/ 38 /M_ 6 = 160 ± 40. If we assume that the covering fraction 
is approximately the volume filling factor, then the average mass-loss rate of the wind, 
M = (1 - f)M + fZM, then leads to L 38 /M-6 = 0.7 ± 0.4. 

We note that the best-fit value for the photospheric optical depth in the blue doublet 
component, ib = 0.62 ± 0.02 is higher than the values obtained for similar stars by fits to 



IUE data (Groenewegen & Lamers 1989). 



4.3. Alternate Velocity Laws 

As discussed in § |3.2j , the change in the maximum velocity of the P Cygni line absorption 
over orbital phase allows us to infer the radial velocity law in the wind, assuming only that 
the X-ray source ionizes the entire wind outside of the shadow of the normal star. In § [4.1| . 
we modelled the P Cygni line variations allowing the X-ray photoionization of the wind to 
vary but assuming a standard (3 = 1 radial velocity law. If we allow a slowly accelerating 
wind with (3 = 1.6 we can still obtain a good match with the observed line profiles, with 
X 2 = 4.91. In this case we find Voo = 1560 ± 35 km s -1 , L 38 /M_ 6 = 0.23 Allowing an 
inhomogeneous wind improves the fit to \ 2 = 4.59. We conclude that while the line profiles 
can be adequately fit with a (3 = 1 velocity law, a more slowly accelerating wind, inferred 
from simpler models, is still consistent with the data. Velocity laws with (3 > 2 (as in some 
of the fits in §|3.2|) generally predict steeply peaked red-shifted emission (Castor & Lamers 



1979); however this is difficult to rule out from the data, as photospheric absorption is 
probably present at low velocities. 



4.4. X-ray Shadow of the Accretion Disk and Gas Stream on the Wind 

The accretion disk or gas stream from the LI point could presumably shadow regions 
of the stellar wind from X-ray illumination. We have made models to find signatures of 
these effects in the P Cygni line profiles. 

LMC X-4 shows a ~ 30 day variation in its X-ray flux (Lang et al. 1981), usually 
attributed to shadowing by an accretion disk that we observe at varying orientations. 
To allow the disk to shadow regions of the wind from the X-rays, we make the simple 
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assumption that it is flat but inclined to the orbital plane. Howarth & Wilson (1983) used a 
similar geometry to model the optical variability of the Hercules X-l system, which shows a 
similar 35 day X-ray period. We allow as free parameters the angular half-thickness of the 
disk, its tilt from the orbital plane, and its "precessional" phase within the 30-day cycle. 
Fixing the 30-day phase to be the at the X-ray maximum (as suggested by observations 
with the All Sky Monitor on the Rossi X-ray Timing Explorer, Paper I), we perform a \ 2 
fit and find a disk half-thickness of 6.9 ± 0.3° and a disk tilt of 31 ± 2° from the orbital 
plane, for a reduced \ 2 — 4.62. 

To simulate the shadow caused by the gas stream from the LI point, we have assumed 
that the stream emerges at an angle of 22.5° from the line to the neutron star (Lubow & 
Shu 1975). The stream blocks the X-rays from a uniform half-thickness p in elevation from 
the orbital plane (as seen from the X-ray source) and at angles < rj from the line to the LI 
point along the orbital plane; p and r\ are free parameters of the fit. The fit gives x 2 = 4.04, 
with p = 33 ± 2° and rj = 90 ± 3°. The values for the other free parameters similar to those 



found in 54.1 



The reason that the gas stream shadow improves the P Cygni profile fits is that at 
0.1 — 0.4, the stream shadows the receding wind, increasing the red-shifted emission, 



which our basic model of § |4.1| underpredicts. Allowing an inhomogeneous wind as in j fP 
improves the fit for the same reason. It is possible that both effects are combined, and 
future observations will be needed to disentagle them. 



4.5. C IV Lines 



The C IV AA 1548, 1550 lines did not show as much variability with the GHRS as the 
N V lines, although we did not observe the lines before = 0.16. Most of the absorption in 
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these lines is probably either interstellar or photospheric (Vrtilek et al. 1997). Nevertheless, 
as a check on our models of the N V lines, we also performed fits to the C IV AA 1548, 1550 



lines. We fix the parameters found in § [4.1| , except for T, 7, and Tb, which are free parameters 
defining the optical depth of C IV in the wind and photosphere. We added fixed Gaussian 
interstellar C IV AA 1548, 1550 absorption lines to the profiles. 

The result of our model, shown in Figure 12, does not match the line profiles in detail 
(x 2 = 11.3). It is possible that the fit is poor because there are photospheric lines in the 
vicinity of 1550A due to ions other than C IV. A systemic velocity closer to might bring 
the photospheric lines into closer agreement with the data. In both the observed profile and 
in the model, red-shifted emission is not prominent. The fitted value of 7 = 7.9 ±0.8 implies 
that C IV in the undisturbed wind is confined to x < 2, so that much of the red-shifted 
C IV emission could be occulted by the primary. If we fit the line profiles with models with 
higher values of (3, then 7 is reduced, but our fits never gave 7 < 4.5 for (3 < 2.4 (the upper 



limit found in §3.2) 



We compare our values of 7 with the optical depth law found for C IV AA 1548, 1550 
by Groenewegen & Lamers (1989) by fitting the C IV AA 1548, 1550 P Cygni lines of similar 
stars observed with IUE. They parameterize the optical depth in the stellar wind by 

r~(^oo) ai [l-(^oo) 1//3 ] a2 (17) 

where 02 = 7 for a\ = and (3=1. For HD 36861 (spectral type 08 III) Groenewegen 
& Lamers find a\ = 1.5 ± 0.2 and 0.2 = 3.1 ± 0.5, which is inconsistent with 7 > 3.1. For 
HD 101413, they find a\ = —0.8 ± 0.3 and a<i = 1.6 ± 0.4. These values imply a sharply 
peaked concentration of C IV towards the stellar surface, but do not fit our data as well as 
7^5. 

We have no compelling alternate explanation for our high value of 7. However, the 
C IV AA 1548, 1550 lines in LMC X-4 that we observed were weak and were dominated by 
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photospheric absorption. Further observations at phases when C IV AA 1548, 1550 is more 
prominent are needed to test our conclusion that this ion is concentrated near the stellar 
surface. 



4.6. Fits using the SEI method 

It has been shown (Lamers et al. 1987) that inaccuracies resulting from the escape 
probability method can be reduced by computing the radiative transfer integral exactly 
along the line of sight while continuing to use the source function given by the Sobolev 
approximation. This is called the SEI (Sobolev with Exact Integration) method. The 
method can be used to simulate the effects of local turbulence within the wind, and using 
the method of Olson (1982), the effects of overlapping doublets can be computed precisely. 

We have implemented a program that uses the SEI method to predict the P Cygni 
line profiles from a wind ionized by an embedded X-ray source. The details of the wind 



ionization are identical to those given by Equations [11] through [Tfj. We find that in our 
analysis of LMC X-4, the approximation we have made for the doublet overlap gives results 
very close to those given by Olson's method and the SEI method. This may result from the 
small amount of the doublet overlap in the N VAA 1238.8, 1242.8 lines, given a separation of 
960 km s _1 between the doublet components and a wind terminal velocity of 1350 km s _1 . 
Using a turbulent velocity v t that is constant throughout the wind, we find that the line 
profiles are best fit with v t < 200 km s _1 . 

The SEI method also confirms the results of the escape probability model for 



the C IV AA 1548, 1550 lines, as reported in § [4.5| . Although the separation of the 

C IV AA 1548, 1550 lines is only 500 km s _1 , we suggest that the doublet overlap was not 

severe during our observations because much of the wind at high velocity had already been 



-28- 



ionized, and C IV is concentrated near the primary's surface, at low wind velocities. 

5. Discussion 

We have demonstrated a new method for inferring the radial velocity profile in a stellar 
wind. The wind radial velocity law for LMC X-4 can be fit with (3 ~ 1.4 — 1.6, and (3 < 2.5. 
This differs from the expected value of (3 — 0.8, but the wind in the LMC X-4 system is 
likely to deviate from that of an isolated O star. Assumptions of the method, that the wind 
is spherically symmetric, for example, will need to be revised in light of further observations 
that cover the entire orbital period. 

Unfortunately, we did not observe the UV spectrum during the X-ray eclipse when 
the X-ray ionization of the blueshifted absorption should be minimal, and so our inference 
of the wind terminal velocity is model-dependent. Observations with STIS and the Far 
Ultraviolet Spectroscopic Explorer (FUSE) during eclipse may distinguish between our 
alternate models. The wind terminal velocity we have inferred from our (3 = 1 escape 
probability method, i>oo = 1350 ± 35 km s -1 , is lower than the terminal velocities measured 
for similar stars by Lamers, Snow, & Lindholm (1995). Their stars # 33,37,38,39, and 44 
are all near the 35,000 K temperature of LMC X-4, but have = 1500 - 2200 km s -1 . 
However, stars in the LMC are known to have terminal velocities ~ 20% lower than their 
galactic counterparts (Garmany & Conti 1985). 

We note from the IUE data shown in Figure 4 that the P Cygni line absorption at 
orbital phases > 0.5 is greater than that at < 0.5. X-ray absorption dips are frequently 
seen at ~ 0.8, possibly indicative of dense gas in a trailing gas stream or accretion wake. 
However, it may be difficult for models to reproduce an accretion wake that occults a large 
fraction of the primary at orbital phases as late as 0.9. Another possible explanation is that 
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a photoionization wake is present (Fransson & Fabian 1980, Blondin et al. 1991). Such 
a wake results when the expanding wind that has not been exposed to X-rays encounters 
slower photoionized gas. While a strong photoionization wake is expected for q ^ 1, 
the presence of Roche lobe overflow and a photoionization wake have been found to be 
mutually exclusive in some simulations (Blondin et al. 1991). A further possibility is that 
at 0~ 0.6 — 0.9 the wind in the cylindar subtending the primary's disk is shielded from 
ionization by the trailing gas stream. A simple prediction of this scenario is that there 
should be more flux in the red-shifted emission at <ft = 0.25 than at </> = 0.75. If the primary 
star does indeed rotate at half the corotation velocity (Hutchings et al. 1978) then the 
gas stream may trail the neutron star by more than the 22.5° predicted by Lubow & Shu 
(1975). 

We have fit the changing P Cygni line profiles using an escape probability method, 
assuming a radial, spherically symmetric stellar wind. However, we note that the 
Hatchett-McCray effect is most sensitive to the region between the q = 1 surface and the 
conical shadow of the primary. Outside of this region, the wind should be entirely ionized 
for likely values of M and L x . The region of the wind that our method is sensitive to is 
illuminated by X-rays, and its dynamics may be affected as a result (Blondin et al. 1990). 
Structure associated with the distortion of the primary, such as the gas stream and a 
wind-compressed disk (Bjorkman 1994), may be also be important in this region. 

We have inferred L 38 /M_ 6 = 0.26±0.01 from our escape probability P Cygni line models. 
Our simultaneous ASCA observations found a 2-10 keV flux of 2.9 x 10~ 10 erg s _1 cm -2 , 
which given a distance to the LMC of 50 kpc implies L38 = 0.9 and M_6 ~ 3 (our 



inhomogeneous wind model of § |4.2| would imply M_6 ~ 1). There is some uncertainty in 
applying L 38 = 0.9, as the isotropic X-ray spectrum and luminosity may differ from that 
observed by ASCA. Nevertheless, the derived M is an order of magnitude larger than 
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that reported by Woo et al. (1995), who used the strength of scattered X-rays to find 
M/voo = 10~ 10 M Q yr _1 (km s -1 ) -1 , implying for our best-fit value of = 1350 ± 35 
km s" 1 that = 0.14. Our result for M is on the order of the single-scattering limit 
given by the momentum transfer from the radiation pressure of the primary; M < L opt /fooC, 
where the optical luminosity of the primary L opt ~ 5 x 10 38 erg s -1 (Heemskerk & van 
Paradijs 1989) implies that M < 2 x 10 -6 M yr _1 . However, if the wind is accelerated by 
scattering of overlapping lines, one can have M = 2 — 5L opt /vooC (Friend & Castor 1982). 

Combining our value of M_ 6 3 with our measured wind velocity near the orbit of 
the neutron star (400 km s _1 ) and the orbital velocity of the neutron star (440 km s" 1 ), 
we find that gravitational capture of the stellar wind (Bondi & Hoyle 1944) could power 
a significant portion of the observed X-rays (L 38 ps 0.3). However, we can not rule out 
that the low velocity and high wind density we have measured are attributable not to an 
isotropic wind, but to other gas in the system, such as the gas stream or a wind-compressed 
disk about the primary. 
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Fig. 1. — A comparison of the LMC X-4 N V line profiles with those of Vela X-l and 
4U1700-37. 

Fig. 2. — Comparison of the N V lines in LMC X-4 observed at = 0.42 with (a) the 
same region in the spectrum of /xCol (dot-dashed; the flux and systemic velocity have been 
normalized to that of LMC X-4), and (b) Gaussian lines at the systemic velocity of the LMC 
and with the same width as optical lines attributed to the photosphere of the primary. The 
vertical lines show the rest wavelenths of the N V lines. 

Fig. 3. — The geometry of the Hatchett-McCray q parameter. Surfaces with constant values 
of q bound regions of the wind at different ionization levels. The primary star is represented 
as a circle centered at the origin, and the X-ray source is marked with an 'X'. 

Fig. 4. — Variations in the equivalent width of a) N V, b) C IV, and c) Si IV. The "*" signs 
mark the equivalent widths from IUE observations reported by van der Klis et al. (1982), 
while the "+" signs are based on the data reported in Paper I. The error bars show the la 
error range. Also shown are theoretical models in which the ion is removed from regions of 
the wind with various values of the Hatchett-McCray parameter q. 

Fig. 5. — The geometry used by a simple analytic method to infer the radial wind velocity 
from P Cygni line variations. The viewer is towards the bottom of the page. The circle 
represents the surface of the primary star and the point marked "NS" represents the neutron 
star. V{4>) marks the highest observed velocity of absorption. In (a) we show the situation 
at <fi = 0.21 and in (b) we show = 0.32. 

Fig. 6. — An example of the measurement of V((f>), the maximum velocity of absorption at 
a given orbital phase. We fit the edge of the observed line profile (the dotted line shows 
the smoothed spectrum) to a horizontal line segment for blue-shifts V > V(<f>) and to a line 
segment with slope as a free parameter for V < V(<j)). The vertical line marks the measured 
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value of V(4>) for the 1242. 8A doublet component, determined by where the continuum 
reaches the same level as for the horizontal line segment. 

Fig. 7. — The wind expansion velocity versus radial distance, as inferred from the P Cygni 
line profiles. The "*" signs mark points determined from the absorption profile. For points 
marked with "+" signs, a model of the blue-shifted emission as the red-shifted emission 
reflected about v = 150 km s" 1 has been subtracted from the P Cygni profiles. For points 
marked with diamonds, the red-shifted emission was reflected about v — km s _1 . Points 
marked with "x" were determined from differences of spectra corrected for blue-shifted 
emission. The solid line shows the best fit for the diamonds for < 0.25 (R > 2). 

Fig. 8. — An illustration of a model that fits the absorption line profiles to determine the wind 
acceleration parameter f3. The conical X-ray shadow of the primary intersects isovelocity 
surfaces at v/v^ = 0.2,0.4. These intersections determine images of the primary's surface 
at Doppler-shifted wavelengths corresponding to v/v^ = 0.2,0.4. The contribution to the 
absorption line profile at v/voo = 0.2,0.4 is determined by the fraction of the primary's 
surface that is not in shadow. 

Fig. 9. — The best-fit model (dashed) to the N V absorption line profiles, allowing the wind 
acceleration parameter (5 and the wind optical depth r as free parameters, and assuming the 
X-rays photoionize the entire wind outside the shadow of the primary. We have subtracted 
a model of the blue-shifted P Cygni emission in order to fit only the blue-shifted P Cygni 
absorption. 

Fig. 10. — The results of the best- fit model of the P Cygni line variations. The right panel 
shows the differences between successive orbital phases, for both the data and model. 

Fig. 11. — The results of the best-fit model of the P Cygni line variations, for the 
inhomogeneous wind model. The right-hand panels show the differences between successive 
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spectra. 



Fig. 12.- 



The results of the best-fit model of the C IV P Cygni line variations. 
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Table 1. 
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* max,0 
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£4 
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(km/s) 




(R*) 


(R*) 


0.095 


0.114 
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1260 


1260 


8.82 


17.57 


17.60 
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0.117 
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8.12 
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1110 
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3.79 


7.45 


7.52 


0.128 


0.141 
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1110 
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3.54 


6.93 


7.00 
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0.200 


490,410 
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3.11 
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450,340 


620,720 
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2.22 
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Table 1 — Continued 
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t / a 
Vobs 
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' max,0 






D 4 






(km/s) 


(km/sj 


(km/sj 


IT) \ 


IT) \ 

(R*) 


IT) \ 


0.293 


0.289 
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430,280 
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1.07 


1.43 


1.75 


0.297 


0.293 


370,190 
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1.06 


1.40 


1.72 
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0.296 


300,260 


490,470 
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1.37 


1.70 


0.304 


0.299 
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1.34 


1.68 
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1.04 


1.32 
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1.59 
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0.316 


340,260 


470,310 
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1.02 


1.21 


1.57 


0.326 


0.319 


350,220 


560,400 


630 


1.01 


1.18 


1.55 


a The 


maximum 


observed 


absorption 


velocity in 


the uncorrected 


N V 



line profile; where two numbers are given these are the blue and red 
doublet components, respectively. All velocities are relative to the LMC (at 
+280 km s" 1 ) 

b The maximum observed absorption velocity in the line profile corrected for 

blue-shifted emission 

c The maximum radial wind velocity based on the observed line-of-sight 

velocity 
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X 2 ^oo P V V subtract 

(km s _1 ) (km s _1 ) 

2.79 1800 ±200 2.3 ± 0.6 200 ± 200 NA 

0.89 1320 ±40 1.39 ±0.14 -40 ± 50 

1.94 1400 ± 100 2.4 ±0.8 -140 ± 130 150 



Table 2: Analytic fits to the wind radial-velocity law 
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Parameter 


Meaning 


Value 
(basic model) 


Value 

(inhomogeneous wind) 


T 


Background wind optical depth maximum 


1.13±0.01 


1.00±0.03,700±200 


7 


Background wind optical depth exponent 


0.77±0.04 


2.1±0.1 




Wind terminal velocity 


1350±35 


1280±35 




Velocity relative LMC 


80±35 


15±35 


Laa/dotM-e 


X-ray luminosity/Mass-loss rate 


0.26 ±0.01 


0.7±0.4 


Tb 


Optical depth of photospheric line (blue) 


0.36 ±0.01 


0.62 ±0.02 


xl, dof 


Goodness of fit, degrees of freedom 


5.7, 690 


4.0, 688 



Table 3: Best-fit parameters for the numeric model of the N V P Cygni profile variations 
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Fig. 10.- 
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Fig. 11.- 
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Fig. 12.— 
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